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Abstract 

Pair creation of strings in time-dependent backgrounds is studied from an effective field the- 
ory viewpoint, and some possible cosmological applications are discussed. Simple estimates 
suggest that excited strings may have played a significant role in preheating, if the string 
tension as measured in four- dimensional Einstein frame falls a couple of orders of magnitude 
below the four- dimensional Planck scale. 
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1 Introduction 



In curing the non-renormalizability of gravity, string theory introduces a large number of 
heavy states: the excited states of a string. The number of these states increases roughly 
as an exponential of their energy, in contrast to field theories, where the number of states 
rises as a power of the energy. This exponential increase is often described as a Hagedorn 
density of states. The extra states in string theory are conventionally assumed to lie near 
the Planck scale. Our grasp of the non-perturbative dynamics of string theory is far from 
complete, but one well-motivated conjecture is that interactions vastly reduce the degrees 
of freedom of the theory, so that the number of available states scales "holographically" as 
the area of a system rather than its volume. Disentangling this "more is less" paradox, 
elucidating the true degrees of freedom of string theory in a non-perturbative regime, and 
reconciling the renormalization group with holography may all be necessary steps before we 
can give a fully satisfactory string theoretic account of the very early universe. We are a 
long way from achieving this, but certainly there is ample reason to believe that theories 
with extended objects are needed to properly formulate quantum gravity, and that such 
theories often have a Hagedorn density of states over an energy range encompassing a great 
many states. Armed with no more than this, we would like to inquire what the possible 
consequences are for cosmology — in particular, for situations where the massive states may 
be produced through the usual ambiguity of the vacuum state in backgrounds with time 
evolution. We will start with an overview of pair production and the steepest descent method 
for estimating occupation numbers. Then we will move on to a possible application to the 
theory of preheating. 

In regimes of parameter space where a spacetime description gives a good approximation 
of string dynamics, the on-shell constraint for a given string state boils down to a second 
order ordinary differential equation in time: 



where x(t) is the wave function for the state in question. This is also the equation for a mode 
of a scalar field, possibly rescaled by some time-dependent factor to eliminate a x term. The 
function u(t) must be determined for any given background and string state. Quantizing 
strings in a general spacetime background is fraught with difficulties because the spacetime 
equations of motion are already encoded in the conditions for conformal invariance on the 
worldsheet, and we seem to be missing some aspects of string dynamics which give rise to the 
weakly curved spacetime we observe. If we can quantize strings in a particular background, 
then oj(t) will be determined, and it will include contributions from the momentum of the 
string and its excitation state. In the following treatment, we will simply assume that 
the spectrum is known, and that it includes the familiar Hagedorn density of states. We 
are interested in establishing conditions under which the total string production is finite 
(without invoking a cutoff on the Hagedorn spectrum), and in estimating the total rate of 
string production when it is dominated by highly excited strings. For this reason, we take 



x + u (t) 2 X = 0, 
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oj{t) 2 to be large. 

To speak meaningfully about string pair production, it's necessary at least to have an 
asymptotic "out" region where u> is slowly varying, in the sense that time derivatives of ui are 
much smaller than the power of uj with the same dimension. When there is such a region, 
one can compare the infinite order out adiabatic vacuum to the actual quantum state to 
determine occupation numbers for a given string state. What the actual quantum state is 
can be subtle, but if there is an asymptotic "in" region where again u is slowly varying, one 
can follow the standard approach of letting the actual quantum state correspond to the in 
vacuum. Let us consider this optimal situation first and further assume that t runs from 
— oo to +00. 

Two widespread analytical techniques for extracting approximate pair production rates 
from (JTJ) are steepest descent contours, applicable when particle occupation numbers are 
small, and parametric resonance, applicable when oo(t) 2 is an oscillatory function of time. 
We will be mainly interested in the former and will develop it in section |21 We will then 
give an application to the theory of preheating in section where conventionally one needs 
parametric resonance to understand the physics — but as we shall see, the steepest descent 
method is still suitable for discussing production of excited string states. The conclusion 
will be that there is a plausible regime of parameters where strings of some type played 
a significant roll in preheating, but that over-production of super-heavy dark matter is a 
potential problem. 

After this work was complete, I learned of ISO], 1 which overlaps significantly with the 
methodology developed in section |21 



2 The steepest descent method 

The steepest descent method is based on early work by pQ and has been developed by various 
authors. Our treatment will to an extent parallel that of [2] . The key assumption is that the 
occupation number \f3\ 2 for a given mode is always much less than 1. Here, (3 is a Bogliubov 
coefficient for comparing the in and out vacua. 2 Setting 

= „-i f duw(u) + e i f duu{u) ^ 

l2uj{t) J2u){t) 



with the requirement |a(t)| 2 — |/3(t)| 2 = 1, one may recast the equation Q as 

ait) = —e 2i P du »(Wf3{t) (3{t) = — e - 2i P d ™m a {t) . (3) 



*I thank E. Martinec for bringing this paper to my attention. 

2 In the standard analogy to one-dimensional scattering, where time is mapped to position and Q is 
regarded as a time- independent Schrodinger equation with potential — w 2 (i), ft is roughly the reflection 
amplitude. 
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Using (3{t) <C 1 and a(t) « 1, one quickly arrives at the general formula for j3 = (3(oo): 







J dt ^- exp (^-2i J duu(u)\ . (4) 



Assuming spacetime is weakly curved in the out region, that \g tt \ — > 1 there, and that uo 
approaches some constant Woo for any given string state, the density of strings states rises 
roughly as e^^, where Th is the Hagedorn temperature, 3 

Tn = 1 (5) 



2irya'c±/6 

for closed strings, where the string tension is r = 1/ (2-ira') and c± is the central charge of the 
transverse degrees of freedom (c± = 12 for the type II superstring) , assumed here to be the 
same in holomorphic and anti-holomorphic sectors. This exponential behavior is modified 
by a power law that depends on the string theory in question as well as the dimensionality 
of non-compact spacetime. Assuming |/3| 2 can be approximated by some function of u^, the 
total number of strings produced may be very roughly estimated to be 

iVtot ~ / ^oo \f3\ 2 e^' T « . (6) 



In s-brane decay (with c^oo replaced by c^o, as explained in the footnote), the exponential 
growth of states is exactly canceled by exponential suppression of (3 at large u, leaving a 
power law behavior that may or may not converge, depending on the dimension 13 IS] - 
Our estimates will focus on the exponential behavior in more generic circumstances. 

A well-known formal tool for exploring the high-energy properties of (0) is to rescale 
t — > t/a, where a is a dimensionless constant. A naive expectation, which will turn out to be 
close enough to the truth in some interesting cases, is that the on-shell condition for highly 
excited string states is given by the rescaled equation, 

X + a 2 uj(t) 2 x = J (7) 

with a 2 set equal to the excitation level N, and with u> remaining nearly the same for different 
string states. More precisely, the assumption is that without any rescaling of time, highly 
excited string states have uo pa auo with a = \fN and u nearly independent of the string 
state. Then (J7J) is correct, with uo — > u, though for a reason orthogonal to time rescaling. 
Granting such a setup, the total number of strings produced is 

/oo 
da\[3{a)\ 2 e^' TH 

r°° uj ( r l _ \ ^ 
/3(a) ~ J dt — exp\—2icr J duuo(u) \ . 



3 Unusual circumstances might invalidate this description of the Hagedorn density, for instance the bound- 
ary operator in [31 01 E] that gives a given open string state a mass that grows exponentially with time. In 
this particular circumstance, w^, in the discussion above could be replaced by u> = too evaluated at the 
time-symmetric point of the full s-branc solution. 
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In (JHJ), oOoo/Th is a fixed number, independent of the string state. Certainly (JHJ) has been 
arrived at through a series of assumptions that are far from self-evident. However it has 
some nice consequences that we believe are more general. First of all, oj(t) 2 must infinitely 
differentiable on the real axis in order to avoid producing an infinite number of strings. The 
same arguments used in [7] to show that occupation numbers in the adiabatic vacuum of 
order A scale as kj - ^" 1-1 ) can be adapted to show that (3 scales as a~ n for large a when 
there is a discontinuity in the n-th derivative of u)(t) 2 . Probably it is also necessary for uj{t) 2 
to be analytic for real t — we will have more to say about this point later. It is intriguing 
that analyticity is also characteristic of the statistical mechanics of systems with a finite 
number of degrees of freedom, hinting once again that the number of degrees of freedom for 
gravitationally coupled strings is finite. 

To evaluate the expression for (5 in (J1J), a technique based on contour integration and 
steepest descent was developed in At least in simple circumstances, the singularities 
of the outer integrand are poles and branch points, occurring in the complex plane where 
u)(t) 2 = or oo. These singularities are distributed symmetrically on either side of the real 
axis because w(t) 2 is real- valued for real arguments. Singularities arising from simple zeroes 
of w(t) 2 were treated in detail in j2], and through a steepest method the following estimate 
was obtained: . ^ 

P » y exp (-2i J T dtu(tf) exp (~2i jf * dtu(t)j , (9) 

where = r — i/j, is the location of a zero of uj(t) 2 in the lower half plane, and r and \x are 
real. If there are several zeroes, one gets a sum of terms of the type appearing on the right 
hand side of (JHJ). Assuming that one zero dominates, we may roughly estimate 

l/?| 2 ~ flVe-^W. (10) 



The factor of ir in the exponent arises from approximating Reu(t) along the line between r 
and t* by an elliptical arc which passes through zero at t* and has its apex at r. A better 
estimate may be obtained if more information is available for the function u(t) 2 . With the 
estimate (jTUjl in hand, we can return to JHJ) and obtain 



tot 



/oo 



Evidently, this converges provided the exponent is negative. If we wished to compute the 
total energy of the strings created, it would only alter the power law prefactor in the integrand 
of (jllj) . and the criterion for convergence would still be that the exponent is negative. 

We derived (JHJ on the understanding that \g a \ — > 1 as t —> oo. If this is not so, but 
there is still an appropriate asymptotic out region, then we need only replace uj^ by 
(defined as a limit). Then the integral in (Jllj) would converge if 

°<T H . (12) 



rtt ,T, 
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The left hand side acts in some rough sense like a temperature for the background in question. 
Similar results can be established for singularities of uj(t) 2 in the complex t plane: the result 
is some coefficient other than ir in (|12|). We will encounter such an alteration of (J12)) at the 
end of this section. 

The type of result expressed in (fT2"j) provides a good intuitive argument, albeit slightly 
circular, for why uj{t) 2 should be analytic on the real line: zeroes and singularities of Co(t) 2 
in the complex plane have to be far enough away from the real axis for (J12|) to be satisfied. 
So the radius of convergence of u)(t) 2 is finite everywhere on the real line. 

As an interesting class of examples, consider a k = FRW cosmology: 

ds 2 = a{r]) 2 {-drf + dx 2 ) . (13) 
The modes of a conformally coupled scalar with mass m satisfy (JTJ with t replaced by r\ and 



UJ 
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k 2 + m 2 a 2 . 4 Assuming that the metric in IjlHj) is the string frame metric, and that the 
curvature is sub-stringy, the mass spectrum is approximately given by m 2 = N/a', where N 
is the level. 5 For highly excited string states, neglecting k is a good approximation, though 
not a uniform one if a(r)) becomes arbitrarily small in the far past. Let us assume that a{rf) 
approaches nonzero constant values, a_oo or aoo, as rj — > ±oo. The scaling analysis above 
is appropriate, with u{rj) = a(rj) / ' \fa' . Assuming that a single zero 77* = r — ifi of a{rf) 2 
dominates, the condition (|T2|) for there to be a finite number of strings created is 

< T H . (14) 



Specializing further, let us consider one of the classic exactly solvable problems: 6 

a(rj) 2 = A + B ttmh(pr]) , (15) 
The exact result leads to exponentially suppressed particle production for large masses: 

W = smh 2 ^ -a; )/2p] ^ ^_ /p = _ M _„ lp 
sinh^oo/p) sinh(7ru;„ 00 /p) 



where U)(rf) = yk 2 + a(rj) 2 m 2 and u; TOO are the limits of uj{rj) in the far past and future, and 
the approximate equality holds good in the limit where uj±oo 3> p. The zero of a{rf) 2 in the 
lower half plane closest to the real axis is 77* = — ^ + ^ log so from (fTUj) we obtain 



7T 



\P\* w e -a^™-c°/P w here c 1 = -^1 + B/A. (17) 



4 Diffcrcnt choices of the parameter £ in the term £x 2 i? that controls the coupling to the Ricci scalar result 
in finite shifts in uj 2 . As long as £ does not grow too quickly as one goes to more highly excited states, it 
should not affect the analysis at the level we arc working at. 

5 Neither the zero point nor the normalization of TV quite agrees with the conventional definition of 
excitation level in type II string theory. This minor discrepancy is of no consequence as long as we correctly 
keep track of the normalization of Tp. 

6 This treatment is similar to the one in 0. 
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Because < B/A < 1, we have ir/A < c\ < 7r/v8, indicating fairly good agreement with 
the exact result (perfect agreement would be C\ = 1). A criterion of the form (J12)) or (1141) 
emerges immediately from estimating the total number of strings produced from (|11|) and 
(fToT) . only with a factor of c\ multiplying the left hand side. Thus we conclude from this 
example that the steepest descent method gives a reasonable approximation of the criterion 
for finiteness of the number of strings produced. 

Let us conclude this section with a background that includes de Sitter space in a certain 
limit, but has well-defined in and out regions away from this limit. The geometry is defined 
by (HHJ) with 

where A, B, and p are constants. If A > B > 0, this background is qualitatively similar to 
(115)1 . but in the limit where A = B > 0, the metric interpolates between the k = patch of 
dSi in the far past and flat space in the far future. To be precise, for A = B and 77 <C —p, 
we have 

L 2 \~A 
aM 2 « -^(-dt 2 + dx 2 ) where L = J— . (19) 
77^ V 2pr 

Although the background defined by ()TH|) and (J18)) is not motivated by string theory con- 
siderations, it might be interesting as a simple "regulation" of de Sitter space: one can for 
example define an S-matrix for A > B > and then take the limit A — > B to investigate 
aspects of quantum field theory in de Sitter space. 

The singularity of a(r]) 2 in the lower half plane closest to the real axis is a branch cut 
starting at 77* = — if p. (There is also a zero at 77 = — - p v ^ r g f , but this is always further 
from the real axis, and it disappears altogether in the A = B limit). The condition (fPI)) 
becomes T s < T H , where 

Ts = -^j-, = -7= • (20) 
ixpa(r) ny/A 

As remarked earlier, T s is an effective temperature for the spacetime. Interestingly, the 
estimate ()20)) does not depend on B at all — though, clearly, B must control a prefactor on 
the total number of strings produced that vanishes as B — >• A. In this limit, T s as estimated 
in (j2IIJ) differs by v2 from the temperature of the de Sitter horizon in the far past, which is 
Tds — l/(2nL) with L as defined in (fT§|) . The discrepancy is entirely due to the crudeness 
of estimating the integral in the second factor of using an elliptical arc passing through 
zero: evaluating that integral exactly for A = B changes the estimate (J2"U|) so that T s = Tas 
exactly. 

It is tempting to interpret the divergence in (jUJ) that arises when (|10)1 is violated as dual 
to the development of an open string tachyon stretched in the complex plane between the 
singularity point t* and its complex conjugate Such an interpretation has been offered in 
time-dependent backgrounds resulting from well-defined stringy constructions jS], and the 
idea of D-branes in imaginary time has been further developed in [S]. 
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In general, one should consider the possibility that the singularities of u(t) (or, more 
precisely, of the integrand in (jSJ)) in the complex t plane are not pointlike or even branch 
cuts, but cover finite regions of the plane. This could happen, for instance, if we replaced 
the a(rj) 2 in (fTo'j) by a continuous sum over different values of p. A criterion like (fT2|) should 
still apply, where /i is roughly the closest approach of the singular region to the real axis. 

3 Strings and preheating 

A widely studied application of particle creation in cosmology is the theory of preheating (see 
for example jHlEH^B)' whereby coherent fluctuations of the inflaton around its minimum 
lead to an oscillating mass term for another bosonic field, \i through a term in the action 
proportional to <p 2 x 2 - The resulting variation in uj(t) 2 can set up a parametric resonance, 
which produces exponentially growing particle occupation numbers for \. The variation of 
oj{t) 2 cannot be expected to be perfectly periodic, since the universe is expanding and there 
may be significant deviations from pure quadratic behavior in the inflaton potential in the 
region where the inflaton oscillates. This limits the amplification that parametric resonance 
can provide. One must then ask under what conditions parametric resonance is "efficient," so 
that an appreciable fraction of the energy stored in the inflaton oscillations will be converted 
to x particles via preheating. Oversimplifying a little, the answer is that the mass of \ must 
vary by a large factor: the maximum mass in each cycle of oscillation must be many times 
greater than the minimum mass. In such a situation, there is a "broad resonance," [TOj, and 
a variety of complicated effects like back-reaction and rescattering become significant. In 
[TT] . more precise analytic criteria were developed to decide whether a resonance is efficient. 
Let us explore how the criteria for efficient extraction of energy from the inflaton oscillations 
might be altered due to the presence of a Hagedorn density of states. 

Occupation numbers of excited string states should be small even if strings do participate 
significantly in preheating. The reason is that resonances become exponentially narrow, and 
the rate of growth of occupation numbers within a resonance becomes exponentially slow, 
as the as one increases the tree-level mass. There can still be a large total number of 
strings produced because of the competing exponential growth of the density of states. To 
distinguish approximately where the competing exponentials have equal "strength," it is 
sufficient to work with |/5| 2 <C 1 for all individual string modes, except perhaps for those 
which are massless at tree level. Our plan of attack, then, is to use the methods developed 
in section |21 to study production of the massive excited states. 

It is very plausible that the effective string tension, as measured in the four- dimensional 
Einstein frame, varies as the inflaton varies. Let us as usual assume a k = FRW cosmology. 
When the string frame metric is weakly curved, the spectrum is given approximately by 
gH r u) 2 = g*trk 2 + m 2 = g^ r k 2 + N/a', where as usual N is the excitation level. The four- 
dimensional Einstein metric is related to the four- dimensional part of the string metric by a 
conformal transformation: ds\ E = e 7</3//Mpi ' 4 <iSg tr where Mpi j4 is the four-dimensional Planck 
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mass, if is a canonically normalized scalar which is some combination of the dilaton and 
the volume of the internal manifold, and 7 is some constant, presumably 0(1) in generic 
compactifications. 7 We may adjust the zero point of (p so that </? = at the minimum around 
which one oscillates during reheating. Evidently, gf E uj 2 = g^k 2 + e lv ' Mvl ' A N/a' . Using the 
usual ansatz ds\ E = —dt 2 + a(t) 2 dx 2 , we end up with 

uj(t) 2 = k 2 a(t) 2 + e^ t)/M ^. (21) 

Now, ip is not necessarily the inflaton: there are many scalars in four dimensions that specify 
the size and shape of the internal manifold. But it seems safe to assume, at least on grounds 
of genericity, that there is some overlap between <p and the inflaton <fi. If <fi oscillates with 
frequency fz between slowly varying extremes (p m i n and <^ max , then we may approximate (f2~Tf 
by 

. . 9 cui + uj 2 u) 2 - u 2 „ 

uj{t) 2 = — t + — t cos nt 

w 2 2 

uo l = k 2 a(t) 2 + e^' M ^- (22) 

a! 

uj 2 , = k 2 a(t) 2 + e™™* /Mp1 -*— . 

a' 

Neglecting derivatives of slowly varying quantities like Q and uj±, the wave- function of a 
given bosonic string mode, suitably scaled by some time-dependent factor, satisfies Mathieu's 
equation, 

Q^ + A-2gcos2^x = 0. (23) 
with z = Vtt/2 and A and q specified by 

A±2 «=Mi) 2 - (24) 

For (A, q) in special regions of the plane, solutions to (|23|) exist whose average value grows 
with z as e vz , where v is the imaginary part of the so-called Floquet exponent. These 
solutions are a manifestation of parametric resonance, but we will not need them for the 
reasons reviewed above. Instead, after the usual "rescaling" uj — > era), we find ourselves in 
the situation described in equations (|7j)- (fT2j) . In (|2*2*jl - (|2*H) . we have neglected the damping of 
the inflaton oscillations due to expansion and to dissipation of energy into the string states. 
Including these effects will make the first oscillation produce more strings than any of the 
subsequent ones. To be more precise: the zeroes of oj{t) 2 = a 2 u){t) 2 as given in are 

t± n =—±-log ^ _ , 25 



7 The expression ds\ E = e lv ^ Mpl ' 4 d,Sg tT is imprecise because the combination of scalars in the exponent 
may not be an eigenvector of the scalar kinetic operator in four-dimensional Einstein frame. But the over- 
simplification will not matter in the subsequent treatment. 
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but we are going to assume that = t~o — t — %\i makes the biggest contribution to string 
production. In order for the formal expression (jSJ) to converge, (|12j) must be satisfied, which 
is to say 



The right hand side is a quantity of order unity. We should regard it as uncertain by at least 
a factor of 2, given that (J22|l is only a rough guess for uj(t) 2 , and the factor of 1/tt is only 
approximate, as the analysis of the example (fTHj) showed. 

If the bound ()26|) is satisfied with room to spare, then the highly excited string states are 
hardly produced at all. Conversely, if it is violated, a huge number of excited string states are 
produced, rendered finite, if need be, by a cutoff on the Hagedorn spectrum. A cutoff at an 
energy A limits the number of strings produced to something like e k l Tn , which is very large if 
A is considerably higher than T#. So it seems likely that the string production is more likely 
to be limited by the available energy in the coherent inflaton oscillations. Thus a violation 
of ()26|) is very likely to mean that excited string states will play a vital role in preheating, in 
the sense that a significant fraction of the energy in the inflaton oscillations goes into excited 
strings. We obtained (|2b[) from a single oscillation of the inflaton. Further oscillations of the 
inflaton are then presumably rapidly damped, and even if a finite number of zeroes of ui{t) 2 
contribute significantly to the total string production, the result will probably only be to 
change a prefactor on the expression for the total number of strings produced: ([26)1 . which 
comes from the exponent of (jllj) . is less likely to be affected. 8 

Neglecting the logarithm in (J26j) as well as various factors of order unity leads to the 
approximate convergence condition 



where r_ is the minimum string tension as measured in four- dimensional Einstein frame. 



The criterion for there to be copious string production due to coherent oscillations of the 
inflaton is that the inequality (|2*7j) should be violated. This happens in a considerably 
broader regime of parameters than the regime that appears in the standard analysis of 
when preheating. In the standard analysis, in order for the resonance to be "efficient," the 

8 We might worry that parametric resonance could alter l|26|) . But the way that occupation numbers 
\P\ 2 significantly larger than i|l(J|) arise in parametric resonance is by having \a\ 2 significantly larger than 
1, so that the approximation a(t) w 1 that led to |0J fails. This is related to having exponential growth 
e vz ss e unt / 2 f certain mode functions. Occupation numbers of order unity for excited string states should 
result in a huge divergence in the total number of strings produced. Exponentially suppressed occupation 
numbers for highly excited strings means that the approximations that led to (|26|) should be valid. 

One can try to test this self-consistent reasoning further by adapting the estimate of to a Hagedorn 
spectrum of bosonic states. The result is dN tot /dt ~ J°° da u{a)e au) / TH , where v{a) is the Floquet exponent 
for a string state at level tr 2 , as determined by the A and q appearing in a rescaled version of (|24|) . This 
expression leads to agreement with <|26ll up to factors of order unity, but it is hard to evaluate those factors. 




(26) 




(27) 



(The point here is that uj- oc ^ffZ while Q is roughly the inflaton mass at its minimum). 
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minimum mass of the additional scalar x must be many times smaller than the inflaton mass, 
so that during many oscillations of the inflaton, the parameters of the Mathieu equation fall 
in the range 2q < A < 2q + y/bq for some constant b of order unity. In contrast, for us, the 
analogous region of efficient preheating is 2q < A < 2xq for some x > 1. Of course, from a 
purely quantum field theoretic standpoint, it may seem contrived for the masses of an entire 
Hagedorn spectrum to oscillate in unison — but that is where the perspective of string theory 
makes natural things that otherwise would seem very special. 

Three potential problems with the idea that strings could have played a significant role 
in preheating are as follows, in increasing order of seriousness. First, the Hagedorn spectrum 
might be cut off at a sufficiently low scale to invalidate the above analysis. Second, the 
fundamental string tension is supposed to be near the Planck scale, much higher than typical 
values for the inflaton mass, so despite the considerations of the previous paragraph, it may 
still seem quite unnatural to have (|2~7j) satisfied. And third, super-heavy stable particles 
could be produced in unacceptable numbers, some of them with fractional charge. Let us 
examine these issues more closely. 

1. Cutoff on the spectrum: Since we believe that the number of degrees of freedom is 
thinned out dramatically by interactions at high energies, we should introduce some 
energy cutoff A and keep only those states with energy less than A. If A is only slightly 
above the string scale 1/a/o 7 , of course all of our arguments collapse, and © is not a 
good way at all to estimate the number of string quanta created. But in fact, there 
is reason to think that A ^> l/\fa': for example, as noted in ^2], the typical size of 
a string is greater than the ten-dimensional Schwarzschild radius associated with its 
total mass until N ~ l/g'f , where N is the excitation level and g s is the string coupling. 
It would be consistent with the spirit of holography to put A at about the energy scale 
of strings at this high excitation level: A ~ l/(v / « 7 fi , f), which is indeed a fairly high 
scale if g s is even moderately small. 9 Since the total energy that can go into the strings 
is finite in preheating, but scales as an exponential of the typical energy of the string 
modes produced, it might be that total energy provides a process-dependent cutoff 
that is often lower than A. 

2. Heaviness of fundamental strings: in conventional scenarios, the square root of the 
string tension is roughly five orders of magnitude larger than the inflaton mass (10 18 GeV 
as compared to 10 13 GeV, roughly speaking). In order to achieve (|2*7j) without making 
the inflaton mass much bigger than it is usually assumed to be, the string coupling 
needs to become very small during the inflaton oscillations — or, more precisely, the 

9 It is argued in |12j that gravitational self-interaction of strings compresses the spectrum and makes 
the density of states rise faster than in the free theory, until at some high scale string states might be in 
one-to-one correspondence with the microstates of black holes. Such a super-exponential increase in the 
number of states one can actually produce in the way we have discussed in this paper is disastrous, since 
almost no background will have a finite total number of quanta produced. Presumably this is where a total 
reorganization of concepts is needed, perhaps with holography playing a central role. 
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four-dimensional Planck scale, which scales as 



Mpi, 4 ~ (VolCFs) 172 ^^'" 2 , (28) 

must attain values much larger than 1 / \foJ . This is not inconceivable, but it does seem 
contrived. It has been pointed out to me ^3] that an inflaton mass during preheating 
substantially larger than 10 13 GeV might be arranged in string theory. This lies outside 
the scope of the current investigation. 

3. Production of super-heavy stable particles: Unless the compactification manifold is 
simply connected (or unless open strings can propagate throughout ten-dimensional 
spacetime), there will be stable states corresponding to strings wound around non- 
contractible cycles. Such states can have fractional electric charge [Hj, and they would 
be very heavy, with mass m ~ £/a', where I is the length of the non-contractible cycle. 
The trouble is that there is a Hagedorn spectrum of excited strings with any given 
winding number, and it is not clear that their production is suppressed strongly enough 
compared to strings with no non-trivial topology which can decay into string modes 
which are massless at tree level. Excited string states with winding number would 
relax to stable, super-heavy particles. The bounds on the density of such particles 
are stringent because they would contribute to dark matter density. Tracing back the 
upper bound on dark matter density today to the epoch of preheating leads to the 
conclusion that super-heavy stable particles should comprise a very small fraction of 
total energy at that time — perhaps on the order 10~ 17 |2j. 

We can estimate how rapidly such particles would be created as follows. The energy 
of states with non-trivial winding is 



to 



2 - u>l + m 2 , (29) 



where uJq is the energy of a string state with no winding, but with the same momenta 
in the non-compact dimensions and at the same excitation level. The density of the 
states with winding should be roughly the same, as a function of lo , as the density of 
states without winding. 10 Thus the calculations leading up to (fTTjl can be adapted to 
an approximate treatment of production of heavy charged particles by replacing the 
density of states 



m 



If production of highly excited string states dominates, then roughly an equal number 
are produced with winding as without, simply because \Jlo 2 — m 2 uj for these highly 
excited strings. Prefactors on the density of states might suppress the production 



10 This is easily seen to be so in the case of a circle compactification when there is winding number but no 
Kaluza-Klein momentum. I doubt the Calabi-Yau case would be much different, but I have not carried out 
the computations explicitly. 
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of strings with winding by a few orders of magnitude — but suppression by a factor 
of 10" 17 seems difficult. This is the relevant factor because we are operating on the 
assumption that excited strings do suck out a significant fraction of the energy of the 
coherent oscillations of the inflaton field. Naturally, if this assumption is lifted, there 
is no pressing problem with overproduction of super- heavy stable particles. 11 

The considerably more stringent limits on fractionally charged particles have not been 
invoked in the discussion of the previous two paragraphs because it's not clear to me 
that fractionally charged particles are a universal aspect of string models with non- 
trivial topology in the extra dimensions. 

In the analysis following (|30[). we have assumed that string production is not somehow 
limited to modes less energetic than £/a'. Intuitively speaking, the point is that the spectrum 
of strings with winding defines the same temperature (despite the presence of a gap) as the 
strings without winding. We have also neglected the annihilation of strings wound in one 
direction with strings wound in the other. Some estimates of the cross-section might be made, 
but annihilation processes don't seem likely to generate a suppression factor anywhere close 
to 1CT 17 . Thus we regard the overproduction of winding states as a major peril if highly 
excited strings are assumed to play any role at all in preheating. 

4 Non-perturbatively constructed strings 

Points 2 and 3 from the end of the previous section may incline us to doubt that excited 
strings can plausibly have anything to do with preheating: on one hand, they're probably 
too heavy to be produced in significant numbers, and on the other hand, if they are, there 
is a tendency to overproduce super-heavy stable particles. 

Strings in four dimensions arising from branes wrapped on cycles of an internal manifold 
offer the possibility of ameliorating the problems described above. Regarding problem 2, the 
key point is that strings coming from wrapped branes become tensionless at special points 
in the moduli space which are a finite distance from generic points and typically occur at 
real codimension two. Regarding problem 3, the helpful feature is that degenerating cycles 
can occur in an isolated nearly singular region of the internal manifold, so that no state with 
non-trivial topology is ever light. We will expand on these points somewhat in the following 
paragraphs and in the process do another simple estimate of string production. 

First, to see that wrapped branes must be included on an equal footing with other states 
in the spectrum, recall for example the duality between the heterotic string on T 3 and M- 
theory on K3, where when cycles of the K3 shrink, wrapped M2-branes become light and 

11 In fact, if we assume that the minimum string tension (as measured in four-dimensional Einstein frame) 
is about an order of magnitude higher than the inflaton mass, then the factor e _7r/iw ' r - 1 suppressing string 
production (c.f. can be arranged to be around 10~ 17 . This might suggest a new twist on the idea of 

string-motivated super-heavy dark matter. 
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provide the W-bosons of enhanced gauge symmetry that can be seen by other means on the 
heterotic side [T5| fTHj) . 

Arbitrarily light strings can arise in various ways [T3 EH 120] • A typical situation is 
for the effective tension in four dimensions to be 

r cS « M\<p\ with M ~ , (31) 

where is a canonically normalized complex scalar field. The precise dependence of M on 
g s may depend on details, but some inverse dependence on g s is to be expected from states 
that descend from string solitons. Let us consider D3-branes wrapped on a shrinking S 2 
within a Calabi-Yau manifold as a definite example. The complex scalar if in this case is 
proportional to $s 2 {Ji + The real part of this integral is the volume of the S 2 , and the 
imaginary part is an axion whose associated instanton is a fundamental string wrapping the 
S 2 . 12 The action includes the standard DBI term, 

S = -r D3 J ^G~±B^+... . (32) 

Integrating this on the S 2 and using td3 ~ l/(a' 2 g s ) gives the Nambu-Goto action for a 
string extended in other directions, with a tension given by (J31|) . The overall normalization 
of the tension definitely depends on more data than we've specified so far (for instance, the 
total volume of the Calabi-Yau manifold), but (pHj) should be right up to factors of order 
unity provided one is not in a peculiar corner of the moduli space (such as a Calabi-Yau 
manifold which is many times larger than the fundamental string scale). 

An S 2 can indeed degenerate at a finite distance in moduli space in a Calabi-Yau com- 
pactification, and in an isolated part of the Calabi-Yau manifold. A famous example is the 
resolved conifold j2lj |22l E3j . It is also possible for degenerate spheres to arise along a two- 
dimensional locus within the Calabi-Yau manifold: this locus might be a Riemann surface of 
non-trivial topology, which resurrects the potential hazard of overproduction of super-heavy 
stable particles. 

At least in certain circumstances, it is clear that non-perturbatively constructed strings 
have the same Hagedorn behavior as fundamental strings. The cleanest evidence comes from 
NS5-branes, where for many coincident branes, the Hagedorn temperature of the fractional 
instanton strings is exactly the Hawking temperature for the near-extremal supergravity 
solution |21]; and from NCOS strings, where the effective string tension of open strings is 
lowered by applying an electric field, and the Hagedorn behavior is simple to understand 
|25j . In both cases, Tjj ~ y/r c g. We will assume that this relation continues to hold, and 
that the arguments that led to (jBJ) still apply. 

12 There are well-understood worldsheet instanton corrections to the metric on the Kahler moduli space 
which we will be neglecting. These singularities are logarithmic and do not change the fact that the point at 
which wrapped D3-branes become tcnsionless is at a finite distance from generic points. They may slightly 
alter (|3TJl . 
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Assume now that the scalar </? has some overlap with the dilaton, and that it is undergoing 
coherent oscillations just after the end of inflation. It is unnatural to suppose that </? passes 
exactly through zero, but it might reasonably be assumed to pass near zero in the complex 
plane. Perhaps naively, we will assume that ip is nearly constant while \<p\ is small. Let ipo 
be the closest approach of ip to the origin. Then the effective string tension varies with time 
like this: 

r cff « ik% + &\ = M^ | 2 + (M*) 2 . ( 33 ) 

Proceeding in the same spirit as in the previous sections, we suppose that modes of the string 
in four dimensions have a dispersion relationship w(t) 2 = k 2 + Nr c g, where N is something 
like the level, and we are not trying to keep track of the expansion of the universe during the 
brief moment when the effective string tension becomes small. The Hagedorn temperature 
is Th ~ y/r e g. The frequencies u(t) (and also the Hagedorn temperature) are slowly varying 
in the far future and the far past, in the sense that derivatives of these quantities are much 
less than the appropriate power of the undifferentiated quantities. So asymptotic in and out 
vacua can be defined, at least in the approximation where we neglect what happens when the 
scalar rolls out of the realm of validity of (|31j). and particle production can be estimated in 
the same way that led to (JTDJ). If we define u(t) = jT e s(t), then the result is that an integral 
of the type found in (JTTj) converges provided a criterion similar to (|26|) holds. Dropping 
various factors of order unity, that convergence criterion is u(r)fi ?Z 1. The solution to 
Tes{t) = in the lower half of the complex plane is t* = r — %\i = — i\ipo\/\ip\. Thus the 
convergence criterion is JM\<po\ 3 /\<p\ 1. We can roughly estimate (p ~ Mpi^minflaton, on 
the assumption that the coherent oscillations of (f are of amplitude Mpi j4 and of frequency 
^mflaton- Recalling that the minimum tension is r c g_ pa M\ip \, we can summarize the 
convergence criterion as follows: 



1 ^ U[r)U = ; = —V ~ TTT7 ~ 10 - „ 

if M<p MMp li4 m inflaton V M pi,4/ (34) 

or equivalently, y/r cS _ ^ 20 Mp1 - 4 ' 

where we have assumed M ~ 10M P14 and m inflaton ~ 10 _5 Mpi j4 . The estimates leading to 
()34|) are inevitably fairly crude as long as we have not specified what the inflaton is and how 
much it overlaps with the scalar tp controlling the effective string tension. The result in (J3*4^ 
is rather different from (|26|) and (|27j). Intuitively, this is because the harmonic dependence 
that led to (|26|) is the smoothest interpolation between a given maximum and minimum in 
a specified time interval, whereas (jSHJ), though still analytic, implies almost a discontinuity 
in the first derivative of u(t) 2 . 

Evidently, the non-perturbatively constructed strings do not have to dip very far below 
the Planck scale in order to violate the convergence criterion (|34[). As before, the exponential 
nature of the divergence then would seem likely to be efficient in sucking energy out of the 
coherent motion of the inflaton. More precisely, it would be efficient in sucking out the energy 
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from the coherent motion of tp, which is probably not the only component of the inflaton. 
This might be modeled in a Hartree approximation by a damping term for the evolution of 
one scalar which acts only near a certain value of that scalar, while other scalars, coupled to 
the first by terms in the potential, are only damped by the usual 3H<p term. 

5 Conclusions 

The convergence criterion on the total number of strings produced, which we have variously 
stated as (fT2"j). and (J3JJ), amounts to the statement that in the appropriate string 

frame, where the string tension is by definition constant, curvatures are substringy — only it 
is crucial that not only the second derivatives of the string metric are bounded by the string 
scale, but all derivatives are bounded in an appropriate way to ensure a uniform radius of 
convergence, of order the string length, for the time-varying energies of string modes. In 
the context of preheating, a violation of this bound on curvatures, either for fundamental 
strings or for strings constructed as wrapped branes, has the simple interpretation that the 
production of highly excited string states efficiently drains the energy of coherent inflaton 
oscillations. This is somewhat in analogy with the standard theory of preheating, but copious 
string production occurs in a larger region of parameters (as measured in the (A, q) plane 
of parameters for Mathieu functions). Unfortunately, it is not particularly clear to me what 
distinctive experimental signatures should be expected if highly excited strings played a 
significant role in preheating. The strings decay into states that are massless at tree level, 
and presumably these can then thermalize. However, if excited strings states are created in 
appreciable numbers, it seems that stable super-heavy particles are likely to be overproduced, 
if they exist in the string spectrum. This problem can be avoided by considering strings in 
four dimensions that arise as branes wrapped on cycles of an internal manifold — cycles which 
shrink to produce a point-like singularity on that manifold. 

The tension of such non-perturbatively constructed strings vanishes at certain points in 
moduli space, and near such points the dependence of the tension on canonically normalized 
scalars may reasonably be assumed to be (|3T|. This dependence leads to a convergence 
criterion (|34j) which is considerably more restrictive than the one derived in (|26|) and (1271) 
assuming harmonic dependence — meaning that the minimum tension doesn't have to be as 
low, assuming the dependence (|31|). in order to violate the convergence criteria and produce 
a large number of excited strings. In fact, as we see in the effective string tension only 
needs to dip a couple of orders of magnitude below the four-dimensional Planck scale in 
order to have copious string production, with otherwise rather standard assumptions about 
the amplitude and frequency of inflaton oscillations just after the end of inflation. 

In summary, it is quite plausible that excited strings (probably of non-perturbative origin) 
play a role in post-inflationary cosmology. 

A much-studied scenario in which light strings appear is ekpyrosis Ej , i n which either 
an M5-brane wrapped on a holomorphic cycle of a Calabi-Yau manifold collides with the E 8 
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boundary of spacetime constructed in [21] , or the two E 8 boundaries collide with one another. 
In the first case, M2-branes stretched between the M5-brane and the boundary behave as 
strings that become light at the moment of collision. In the second case, M2-branes stretched 
between the two E% boundaries become light, and they become precisely the perturbative 
heterotic strings in the limit where the two boundaries coincide. The results of this paper 
bear on ekpyrosis in two ways: 

1. We have noted in section 0] a natural way of avoiding or "regulating" the collision: one 
simply needs to assume a non-zero value for the pseudoscalar superpartner of the real 
scalar controlling the distance between the branes/boundaries. In the case of colliding 
boundaries, this peudoscalar is just the axion arising from the Neveu-Schwarz field 
of the heterotic string with both indices in four dimensions. If this axion exists in the 
spectrum, it is in fact unnatural to suppose that it is zero during ekpyrosis. 

2. If ekpyrosis does occur, presumably in a "regulated" way as explained in the previous 
point, then estimates of the creation of light strings proceed in parallel to our cal- 
culations in section |3J If excited strings are produced during ekpyrosis, the problem 
of overproducing stable super-heavy particles seems likely to recur. Specifically, if an 
M5-brane wrapped on a holomorphic curve hits an E 8 boundary, then the light strings 
can wrap cycles of the holomorphic curve (which is a Riemann surface embedded in the 
Calabi-Yau manifold), unless that curve is topologically an S 2 . And if the two bound- 
aries collide, then the light strings can wrap on the Calabi-Yau, unless the Calabi-Yau 
is simply connected, which would considerably reduce the particle-phenomenological 
appeal of the setup. 

It doesn't help to assume that the strings are light only for a brief instant of time. 
In fact, according to (jHH), larger ip (meaning a briefer collision) causes more string 
production, not less: it's then easier to violate the inequality in (J21J). 

The reason that it is easy to offer the plausible interpretation of draining energy from 
coherent inflaton oscillations when strings are copiously produced just after inflation is that 
a significant fraction of the energy is kinetic (relating to the time-derivative of the inflaton). 
A much thornier question is what happens if a string becomes lighter than the Hubble scale 
during slow-roll inflation, when most of the energy is stored in the scalar potential. I have 
speculated elsewhere j2H] that light strings might limit the rate of inflation to the string 
scale, and that there could be a fixed point of forward time evolution where a string tension 
is vanishing at the same time that the expansion of the universe is slowing to something less 
fast than exponential. An observation that tells against this idea is that string production 
should result in a new source of positive energy density, and unlike the case of preheating, 
there is no obvious way to suck this energy away from the scalar potential energy. Granting 
the Friedman equation H 2 = pj (3Mpj 4 ), it is hard to see what mechanism due to light strings 
will drive p and therefore H to zero. Perhaps this question must wait on some treatment 
more firmly based on the first principles of string theory. 
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Indeed, an objection that could be raised to the entire enterprise of this paper is that 
a' corrections to the spacetime equations of motion become significant when curvatures are 
string scale, so it's not clear that we have adequately addressed classical effects before delv- 
ing into the quantum effect of particle creation. There are at least three reasons to continue 
thinking about string creation in general backgrounds: 1) Non-perturbative strings don't 
have a well-understood connection to spacetime equations of motion. 2) Worldsheet tech- 
niques have not so far provided a particularly large class of well-understood time-dependent 
backgrounds. 3) Because copious string creation occurs, in some heuristic sense, at the 
same order in a 1 as corrections to the classical spacetime equations, string production prob- 
ably does have an important role to play in our understanding of general time-dependent 
backgrounds. 

Let us end with some open questions: 

• Our estimates for string production were fairly crude in places, starting with the neglect 
of power law corrections to the exponential growth in the density of states, and includ- 
ing also the neglect of various factors of order unity that enter into the exponential 
behavior in Is it possible to give more precise estimates? 

• Can one give an explicit string theory construction of time-dependent backgrounds in 
which highly excited string states are produced? This is surely complicated by the fact 
that back- reaction from such strings is likely to be important. 

• Might the CP violation that results from turning on an axion to "get around" a point 
of tensionless strings be of interest in baryogenesis? 

• If highly excited strings are produced when inflation is over or nearly over, is there 
some definite affect on the spectrum of fluctuations? How do decay products of highly 
excited strings thermalize? 

• If we assume that the end of inflation involves oscillations in a complex moduli space 
with tensionless strings on some real codimension two loci, how much isocurvature 
perturbations are likely to be generated? Is the scenario plausibly consistent with 
existing bounds? 

We hope to return to some of these questions in the future. 
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